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I^H ' Abstract 

' Gravitation lensing calculations, which are generally done for light 

, ray, are extended to that for a massive particle. Many interesting results 

were observed. We discuss the scattering cross section along-with many 
consequential quantities here. In particular, the case of Schrwarzschild 
metric was taken as illustration, though the analysis applies to a wide 
^— ^ ' range of cases, such as extended black holes. 

o : 

,—1 ■ 1 Introduction 

m ; 

■ Recently, there has been a lot of interest toward fixing the mass of a neutrinopP |2] , 
^ I which constitute a large part of the cosmic background. It is anticipated, that 
O^' such a neutrino flux can be used to probe the gravitational potential of the scat- 

I I tering object. In conventional relativity text books, a calculation of such a pro- 

&J!)' cess is done with the case of light ray (particle with vanishing rest mass) , which 

^ I is the simplest possible, though very important case of gravitational scattering. 

• '~j ■ Even though, for massive neutrinos moving with large momentum, light-like ap- 

r> I proximation should work with no much deviations from the actual experimental 

■ results, the case we discuss is also of importance, since it is observed that there 
is no proffered range of energies with which cosmic background particles come 
in and hence, one should not a priori be driven by the assumption of getting 
high energetic particles always. 

In the present paper, we discuss the various quantities relevant in experimen- 
tal observations of such a scattering process for a massive particle. In particular 
we define two quantities, namely magnification ratio (/i) and refractive index 
(77). Following standard notations, magnification ratio, in general, can be de- 
fined as the ratio of the apparent area, projected on the source plan, of the 
measured flux to the area of the incoming flux on the same plane. In this work, 
we try to attach a refractive index to the space, following similar ideas from the 
usual optics. Refractive index, is defined as the proper instantaneous velocity of 
the particle at a spatial point to its asymptotic value. It can be anticipated, that 
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since defining above quantities use of only geometric structure of the system is 
made, the quantities are apphcable as soon as metric, and hence geometry of 
the space-time is given. It is not possible to give a general formula for these 
quantities, as can be seen from their metric dependence and/or choice of the 
co-ordinates. We illustrate the idea by deriving these quantities for the specific 
metric we have chosen, Schwarzschild metric. 

Many interesting, new features were observed, compared to the case of light 
scattering. The analysis turns out to be of no much complication, but in- 
structive, in this sense. The main course, of the paper, is divided into various 
sections. In section II we start with deriving scattering results for the case of 
Schwarzschild metric. In section III, we discuss the magnification ratio for this 
case. Refractive index for this metric is studied in section IV. 

2 Scattering Calculations 

The geodesic equation for the present case^ (static, isotropic, single compact 
center) is given by[3j, 



The constants of motion (i.e. E Sz J), bear the same physical meaning as in 
the classical regime, i.e. J is angular momentum and (1 — E)/2 is energy in the 
classical regime. 

The analysis will go on the same line as for the case of light (as per the 
standard textbook calculations), but this time we have more complications. 
Unlike the case for light, we are bound to keep both J and E non-zero. However, 
as mentioned in the introduction, this doesn't make the analysis very compli 
cated and one can easily proceed further. 

Eliminating the parameter 'p' from the above equations and substituting 
Robertson Expressions for A{r) and B{r) (both truncated up to second 
terms, following standard analysisPI), we end up with a relation between the 

co-ordinate and V co-ordinate as 



One could take higher approximations in A{r) without much difficulty. B(r) 
seems to be more troublesome, but it can be seen that, there is no harm in 
including one higher order, which would change our definition of J, in further 
analysis. But, since one expect a very small correction even at the higher 
energies(and even at the speed of light!), there is no need to worry about the 
corrections introduced by higher approximations. 

^here, we mean Schwarzschild metric 




(1) 
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The above integral splits into nice parts, one equal to the classical expression 
and the next as a relativistic correction, 



(4i£_ + ,GMr , t^)* (3) 

^(i-is) + aeH-iJ Jr.... ^(1 _ £) + sew _ ij 

Small algebraic rearrangements give, in the limit r — > oo 



G^M^\ r H)dr GM 



The integral in the first term is just the classical expressiongl. The multiplica- 
tive factor contains the relativistic correction, which is of the order of (G^M^/ J), 
and which decreases as impact parameter increases. This is expected, since far 
away from the source^, space-time is essentially flat. On the other hand, this 
correction blows up at J = limit. It can be seen from the comparisons shown 
in fig. 2, that this quantity cease to be meaningful immediately after it reaches 
approximately the limiting value set by the Schwarzschild radius. One can not 
interpret the result by giving any simple meaning to the results then, since all 
the equation and hence the analysis is invalid in the region given by r < Rsch- 
The comparison between bending of light ray and a massive particle shows that 
there is a clear cut distinction between the two cases (fig. 3) 




Figure 1: fig a. Comparison between classical and relativistic results for a 
massive particle around Sun-like center, fig b. Comparison between Angle of 
bending for light and that of a massive particle 



Having calculated (j)o, it is straightforward to get X: ^^nd hence differential 
cross section da, since the relation^ between differential cross section^ and X) 

^property due to compactness of the source 

•^since, the problem reduces to an equivalent central potential problem 
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is same as that for Coulomb scattering. i.e. 

^- = fc"-^ (5) 
sm (x/2) 

where, k' = GM/{1 — E) and is the solid angle. 

On the other hand, it should be noted that, because* of the possibility of 
multiple images, in general the experimental measurements are more involved. 
Following figure illustrates this point. 




Figure 2: Various components in the flux due to Multiple Images 



The no. of particles arriving the observer per unit time per unit area of cross 
section as O is given by. 



dn = drii = ■ 



dp 
dO 



dn (6) 



To Sin f 

Summation is over each value of p, satisfying 

pi + To sin 6 = [to cos 6 - pi tan(xi /2) ) tan Xi (7) 

since X is a known function of p, the above equation can be solved for allowed 
values of p for given 6 and robs- (Note that the exact expression for the flux will 
contain the components of dui 's in the direction 0, but in the limit robs ~^ cxd 
the corrections are insignificant) 

In the similar fashion one can get the relation between r and t as, 
H , r [l + (2 + 7)^]rfr- 



the situation can be thought as the effect of an attractive potential 
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which again can be separated into two parts, classical and relativistic correction 
as, 



\J (^^ - - 2a(r - r~n) 



(9) 



with now a = (GM)/(1 — £■). Compared to the classical result^ there is a clear- 
cut demarcation, now. Classically, the particle takes less time to travel between 
any two points (symmetric to r„ij„) in presence of the scattering center than it 
would have taken if it were moving on the straight line joining the two points 
when center is now there, whereas relativistically, the particle takes more time. 
This result is quite interesting, since it can be utilized as a check for validity 
of the relativistic domain. In the presence of the field, the particle traveled 
more distance between two spatial points, but with a velocity (of magnitude) 
always greater than v^a i-c. it traveled greater path with a greater velocity. In 
classical case, the additional path introduced due to bending must be such that 
the particle took less time, whereas in relativistic case, there are two differences. 
The path length in this case is more (because of Robertson approximation) and 
the particle has to travel more distance than the classical one^. At the same 
time there is an effect of time dilation, which was absent in the classical case. 



3 Magnification Ratio 

The magnification ratio is defined as the ratio of the projection of the apparent 
area of the flux of the particles in presence of the scattering center to that if it 
were not there. Since, we assume the incoming particle flux to be of uniform 
areal density, it can be shown that this quantity equals, the ratio between the 
apparent area of the flux to area through which the flux was incident. The idea 
is illustrated in fig. 2 

In the case (see fig. 2), when there is cylindrical symmetry (which in most 

^Classically, one gets 



na3 , {r-a) + ^ (r^ - r^.„) - 2a{r - Tmin) 
In 

which is nothing but 

T = t(ri,rmin) + t(r2,rmin) 

Voo 

^Unlike classical case, now the work done on the particle by the field is used to 1. increase 
its speed and 2. increase its inertia 



5 




Figure 3: Apparent path diagram (analogous to ray diagram) for the particle 
in a cylindrically symmetric system. 



cases can be assumed without trouble), it is given by 

h dh 

P dp 

Further, from the geometry of the figure (x = 20o — tt), we have 

h = p + rs tanx(/o) 



(10) 



(11) 



Given metric, x can be expressed in terms of p, the differential ^ and hence /U, 
is well defined (since x is a monotonically varying function of p). 
In our case we get. 



P- 
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A{k'/pf{rjp) ^{k'/pf{rjpf 



(12) 



(1 - {k'/pY) (1 - {k'/pYY 

where, k' = (GM/w^). The relation between differential cross section and 
magnification ratio can be given by, 

2k' 



da = — 

P- 



(l-tan^(x/2)) 



+ rstan (x) 



dr. 



(13) 



The above result is expected, since the differential cross section is a function of 
boundary conditions and possibly of metric only, whereas /i does depend upon 
Tg, the distance between object plane from the scattering center. 
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4 Refractive Index 



We define local refractive index of a spatial point as the ratio of the magnitude of 
the asymptotic velocity of the particle (v^.), to the magnitude of the velocity at 
that point. In the classical case, where the above definition is unambiguous, the 
relativistie case needs more specification, as there are more than one velocities. 
It can be shown that it is easier to work with the proper velocity of the particle. 
Hence, 

rj{x, y, z) = ^ (14) 

From the symmetry of the Schwarzschild metric, it can be guessed that the 
refractive index will not depend upon cj) and 6 co-ordinates. The main concern 
is, whether it is possible to determine the trajectory of the particle completely, 
once the initial conditions are supplied. Because, only then one will be able 
to talk about the equivalence between specifying the metric and specifying the 
refractive index. 

We start with the classical result, 



-f 

Jr„ 



{J/r'^)dr 



J2 



(15) 



J/r- 


f mk/ J 






^2m{E - 





which under proper b.e. leads to the expression, 

J/r + mk/^ 

cos^= — , ... , (16) 

Hence, 

sin^=-5! — ' ' =J— (17) 

Now, if we define ry(r) as, 

sin((/. oo) 
sm((/)(r)) 

which is nothing but the typical Snail's law in optics, then from eq.l4, it can 
be seen that this definition matches with our definition in terms of the local 
velocity. 

For relativistie case, one would go along the same line and it is again easy 
to verify that our definition is equivalent to that obtained by Snell's relation. 
And from our experience, in optics, that once initial angle of incidence is given 
for light ray, using Snell's law one can determine the path completely, one can 
deduce that, once impact parameter of the particle is given, which is equivalent 
to specifying the angle of incidence, one can use the refractive index to determine 
the trajectory. In case of particle, one more information must be provided, 
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namely incident velocity (which is equivalent to specifying momentum and hence 
in case of light ray, specifying its wavelength, which will determine the opacity 
of the material). 

In this case, one can define the refractive index as 



.^^^l/B{r)-J^/r^-E 



5 Conclusions 

From the above treatment, it can be deduced, that, that neutrinos, if found 
to be massive, will follow a trajectory different from that of light ray will, 
corrections depending upon the value of the impact parameter. The time delay 
measurement done accurately enough, can be used to test the validity of the 
energy domain for application of classical results. On the other hand, the later 
analysis shows that, the analogy between gravitational lenses and optical lenses 
can be extended further to include refraction effects. 
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